XMM-Newton observations of the Perseus Cluster I: The 
temperature and surface brightness structure 
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ABSTRACT 

We present preliminary results of the XMM-Newton 50 ksec observation of the Perseus cluster 
which provides an unprecedented view of the central 0.5 Mpc region. The projected gas temper- 
ature declines smoothly by a factor of 2 from a maximum value of ~ 7 keV in the outer regions 
to just above 3 keV at the cluster center. Over this same range, the heavy element abundance 
rises slowly from 0.4 solar to 0.5 solar as the radius decreases from 14' to 5', and then rises to 
a peak of almost 0.7 solar at 1.25' before declining to 0.4 at the center. The global east/west 
asymmetry of the gas temperature and surface brightness distributions, approximately aligned 
with the chain of bright galaxies, suggests an ongoing merger, although the modest degree of 
the observed asymmetry certainly excludes a major merger interpretation. The chain of galaxies 
probably traces the filament along which accretion has started some time ago and is continuing 
at the present time. A cold and dense (low entropy) cluster core like Perseus is probably well 
"protected" against the penetration of the gas of infalling groups and poor clusters whereas in 
non-cooling core clusters like Coma and A1367, infalling subclusters can penetrate deeply into 
the core region. In Perseus, gas associated with infalling groups may be stripped completely 
at the outskirts of the main cluster and only compression waves (shocks) may reach the central 
regions. We argue, and show supporting simulations, that the passage of such a wave(s) can 
qualitatively explain the overall horseshoe shaped appearance of the gas temperature map (the 
hot horseshoe surrounds the colder, low entropy core) as well as other features of the Perseus 
cluster core. These simulations also show that as compression waves traverse the cluster core, 
they can induce oscillatory motion of the cluster gas which can generate multiple sharp "edges" , 
on opposite sides or the central galaxy. Gas motions induced by mergers may be a natural way 
to explain the high frequency of "edges" seen in clusters with cooling cores. 



Subject headings: galaxies: active - galaxies: 
individual: NGC 1275 - X-rays: galaxies 

1. Introduction 

The Perseus cluster of galaxies (Abell 426) is 
the X-ray brightest nearby cluster (with z — 
0.01756, 1' corresponds to ~30 kpc (29.7 kpc) for 
H = 50 kins" 1 Mpc -1 which we use through- 
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out). It was first reported as an extended X-ray 
source from UHURU observations (Forman et al. 
1972) and first images (Gorenstein et al. 1978; 
see also Fabian et al. 1974; Wolff et al. 1974; 
Malina et al. 1978) provided a glimpse of the X- 
ray structure of Perseus. Detailed X-ray images 
were obtained with the Einstein IPC (Branduardi- 
Raymont et al. 1981) and HPJ (Fabian et al. 
1981), the ROSAT PSPC (Schwarz et al. 1992, 
Ettori, Fabian, White 1998), HRI (Bohringer et 
al. 1993) and recently with Chandra (Fabian et 
al. 2000, Schmidt, Fabian & Sanders, 2002). The 
cluster has a prominent X-ray surface brightness 
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peak centered on the active galaxy NGC 1275, 
containing a strong core-dominated radio source 
(Per A, 3C84) surrounded by a lower surface 
brightness halo (e.g. Pedlar et al. 1990, Sijbring 
1993). Many studies have explored the correla- 
tions between X-ray, radio, optical, and ultravi- 
olet emission in the core (see e.g. McNamara, 
O'Connell & Sarazin 1996 and references therein). 
The proximity and brightness of the cluster make 
it a natural target for studies of a relatively re- 
laxed cluster with a cold core. In this contribution, 
we report preliminary results of the XMM-Newton 
50 ksec observations of the Perseus cluster. 

The structure of the paper is as follows: in Sec- 
tion 2 we discuss XMM-Newton observations of 
the cluster, in Section 3 we describe the analy- 
sis of the data and present the results, in Section 
4 we discuss possible implications of the results, 
and in Section 5 we summarize our findings. The 
analysis of the substructure in the very core (less 
than 2') of the Perseus cluster and details of the 
heavy element abundance and distribution are be- 
yond the scope of this paper and will be reported 
in subsequent publications. 




Fig. 1. — Surface brightness contours in the 0.3-5 
keV energy range (corrected for the background 
and telescope vignetting) superposed onto a DSS 
optical image (30' x 30'). The X-ray image has 
been smoothed with an 8" Gaussian. The central 
galaxy is NGC 1275. A chain of bright galaxies is 
visible to the West. 



2. Observations and background subtrac- 
tion 




Energy, keV 



Fig. 2. — The spectrum of the 8' — 13.75' annulus 
around NGC 1275 and the background spectrum. 
Two lower curves show the blank field "quiescent" 
background and the additional "flare" background 
component (power law) assumed in the analy- 
sis. Two upper curves show the source spectrum 
with the quiescent background subtracted (upper 
curve) and the source spectrum with both the 
quiescent and flare background components sub- 
tracted (lower thick curve). The normalization of 
the flare background has been estimated using the 
data above ~11 keV. 

The core of the Perseus cluster was observed 
with XMM-Newton on 2001, Jan. 30 with a 
total integration time of ~53 ksec (observation 
ID #00851101). The pointing was centered on 
NGC 1275 (Fig. 1). We report below on the re- 
sults obtained with the MOS instrument only. Re- 
sults from the RGS and PN observations will be 
reported elsewhere. Calibrated event lists were 
generated using SAS v5.3. As is known the MOS 
background has a steady (quiescent) component 
partly associated with cosmic X-ray background 
and variable component due to charged particles 
(e.g. Lumb et al. 2002). To check the level of the 
variable background component, we used the same 
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tracer as implemented in the SAS task emchain — 
i.e. the count rate of events with energy deposited 
in a single pixel higher than a threshold value 
("REJECTED_BY_GATTI" - in SAS notation). 
Although there are no obvious spikes in the light 
curve of such events, the count rate is several times 
higher than the recommended threshold value for 
the entire observation, indicating a significantly 
enhanced background level. Thus a straightfor- 
ward application of the recommended cleaning 
procedure (using emchain default settings) would 
remove essentially all the data. To use the dataset, 
we had to adopt special procedures to correct for 
the variable background component. As a first 
step, the steady component of the background 
was subtracted using data from blank field obser- 
vations, provided by the XMM SOC site (back- 
ground data files are from 2002, Jan., see Lumb et 
al. 2002). The spectra of all but the central chip 
were accumulated for MOS1 and MOS2. Assum- 
ing that the variable component of the background 
has a power law shape (in counts/channel space) 
with the slope of B{E) = A x E a and using the 
data in the lf-f 2 keV range, we calculated the nor- 
malization of the variable component. The slope a 
of the variable background spectrum is known to 
vary from observation to observation and (rather 
arbitrarily) has been fixed at a = —0.45. We then 
assumed that the intensity of the variable com- 
ponent per unit area is constant across the de- 
tector and corrected all observed spectra for this 
variable component. The sample source and back- 
ground spectra for the 8' — 13.75' annulus is shown 
in Fig. 2. The same procedure applied to the pub- 
licly available datasets with the "quiescent" back- 
ground typically yields a normalization for the ad- 
ditional background component that is at least an 
order of magnitude lower. 

To verify the quality of the background sub- 
traction procedure we used a publicly available 30 
ksec observation (OBS1D 0100240801) performed 
<~14 hours after the end of the Perseus observa- 
tion. In this observation the background level was 
still enchanced compared to the quiescent level, 
although a factor of ^2 lower than in the Perseus 
observation. There are no strong sources in this 
field (apart from the target - a compact source) 
and it is easier to see how well the background is 
corrected in the regions far from the center of the 
field of view. For this data set we repeated exactly 



the same background correction procedure as for 
the Perseus observation. The results are shown 
in Fig. 3. One can see that above ^4 keV the 
predicted extra background component matches 
almost perfectly the accumulated spectrum (cor- 
rected for quiescent background). Below 2 keV 
some residuals are present, but given the bright- 
ness of Perseus at these energies (factor 10 to 100 
brighter than the background) such residuals have 
no significant impact on the results of the spectral 
fits. 

The procedure and functions described above 
are the result of extensive tests with various shapes 
of the additional background and different meth- 
ods to determine its normalization. The results 
reported below were found to be robust to the re- 
maining uncertainties in the background subtrac- 
tion. Further refinement of the background cor- 
rection procedure may change slightly some of the 
numerical values, in particular the maximum gas 
temperature close to the edge of the field of view 
by ~5-10%, but cannot affect any of the principle 
results in the paper. More subtle features, like the 
limits on the possible nonthermal emission compo- 
nents, which are in principle attainable with the 
required statistics are deferred to future work. For 
the same reason we report only on the MOS ob- 
servation, for which we achieved the most robust 
results. 

For the subsequent analysis we used MOS data 
with patterns in the range 0-12 and the recom- 
mended value of the flag (XMM_EA). For the 
spectral analysis we use one of the MOS response 
matrixes provided by the XMM SOC (namely 
"ml_thinlv9ql9t5r5_all_15.rsp") and assume that 
the same response (corrected for energy dependent 
vignetting) is applicable for all regions. 

3. Analysis 

Although the cluster is not perfectly symmet- 
ric (see Fig. 1) it is illustrative to derive the ra- 
dial dependence of major parameters, in partic- 
ular the density, temperature and heavy element 
abundance, assuming spherical symmetry with the 
center at NGC 1275. 

The radial dependence of the surface brightness 
is shown in Fig. 4. The dotted line shows the ex- 
pected profile (taking into account the PSF of the 
telescope) for a point source plus constant sur- 
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Fig. 3. — The same as in the previous figure 
but for the observation performed ~14 hours af- 
ter the end of the Perseus observation. The back- 
ground level was still enchanced compared to the 
quiescent state, although a factor of ^2 lower 
than in the Perseus observation. The upper curve 
shows the quiescent background spectrum. The 
histogram is the spectrum accumulated over the 
8' — 13.75' annulus and corrected for the quiescent 
background. The power law shows the predicted 
"flare" background component, obtained with the 
same procedure as was used for the Perseus obser- 
vation. 

face brightness. The central excess (within 10") 
is therefore consistent with being due to the con- 
tribution of the unresolved, compact source - the 
nucleus of NGC 1275. 

3.1. Spectrum of NGC 1275 

Given the XMM-Newton spatial resolution and 
complicated surface brightness distribution near 
the cluster core, the precise determination of the 
NGC 1275 nuclear spectrum is difficult. We used 
the central 14" (radius) circle and the 14"-28" an- 
nulus as the source and background regions re- 
spectively. The resulting spectrum is shown in 
Fig. 5. It is fairly well fit with an absorbed 
(Nh = 1-2 x 10 21 cm^ 2 ) power law with a photon 
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Fig. 4. — Radial dependence of the surface bright- 
ness in the 0.3-5 keV band. The dotted line shows 
the sum of the telescope Point Spread Function 
(PSF) plus a constant, illustrating that the cen- 
tral excess is consistent with being due to the 
NGC 1275 nucleus contribution. 

index of 1.65 and a narrow line at an energy of 
6.289 keV (with la confidence interval from 6.276 
to 6.313 keV) and an equivalent width of ~165 eV. 
The assumption that the line is due to fluorescent 
emission of neutral iron at 6.4 keV implies a red- 
shift of z = 0.0176 ± 0.003, fully consistent with 
the optically determined redshift of NGC 1275 
z = 0.01756. A significant amount of cold gas 
and dust is known to be associated with the high 
velocity (HV) system (Minkowski 1957, see also 
Conselice, Gallagher, & Wyse 2001), but the ex- 
cellent agreement between the redshifts of the line 
and the galaxy indicates that reprocessing by the 
HV gas does not contribute significantly to the ob- 
served spectrum. The 6.7 keV line which is dom- 
inant in the cluster gas emission is very weak in 
the nuclear spectrum, suggesting that background 
has been subtracted with sufficient accuracy. We 
conclude that the observed spectrum of the nu- 
cleus is consistent with typical AGN spectra. The 
luminosity of the nucleus derived from the best fit 
parameters is of the order of 10 43 ergss _1 in the 
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Fig. 5. — NGC 1275 spectrum determined as the 
difference between the spectra extracted from the 
14" (radius) circle and the annulus from 14" to 
28". The position of the line is consistent with 
the neutral iron fluorescent line at 6.4 keV with 
the rcdshift of 0.0176. 

0.5-8 keV range and is subject to systematic un- 
certainties (at the level of 20%) due to background 
subtraction and PSF correction. 

3.2. Projected parameters 

Spectra were accumulated in a set of annuli cen- 
tered at NGC 1275. The resulting spectra were fit 
(in the 0.5-9 keV range) with a single tempera- 
ture APEC (Smith et al. 2001) model in XSPEC 
(version 11.2.0; Arnaud 1996) with the gas tem- 
perature T e , abundance of heavy elements Z (rel- 
ative to the solar values of Anders and Grevesse 
1989) and low energy photoelectric absorption col- 
umn density Nh as free parameters. The red- 
shift has been fixed at the optically determined 
value for NGC 1275: z = 0.01756 (Strauss et 
al. 1992). The projected radial dependence of 
all parameters is shown in Fig. 6. The derived 
value of Nh ~ 1.2 — 1.3 x 10 21 cm~ 2 is compa- 
rable with (and even somewhat lower than) the 
galactic hydrogen column density in this direction 
Nh ~ 1.5 x 10 21 cm,- 2 (Dickey & Lockman 1990). 
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Fig. 6. — Radial profiles of Nn,T e and heavy el- 
ement abundance calculated for a set of annuli 
centered at NGC 1275. The parameters were ob- 
tained using a single temperature APEC model. 
The first two data points in the innermost region 
are affected by the presence of the bright compact 
source - the nucleus of NGC 1275. 

No intrinsic absorption is required by the data 
and no large variations are observed with radius. 
However we note that, given the uncertainty in 
the background subtraction procedure described 
above, it is possible that we systematically un- 
derestimate the absorbing column density. The 
temperature shows a clear decrease from 7 keV at 
large radii to 3.7 keV near the center of the clus- 
ter, except for the very central region which is con- 
taminated by hard emission from the NGC 1275 
nucleus. The abundance of heavy elements shows 
an opposite trend, decreasing (outside the central 
1' — 2' region) with distance from the center. 

3.3. Asymmetric surface brightness sub- 
structure 

As is clear from Fig. 1 and known from earlier 
Einstein and ROSAT observations (Branduardi- 
Raymont et al. 1981), the surface brightness is 
not symmetric around NGC 1275. In order to fur- 
ther enhance the asymmetric structure, we plot in 



5 



Fig. 7. — Relative deviation of the surface brightness from the azimuthally averaged value. The image shown 
is 30' x 30'. Clearly seen are: the region of excess emission to the East from the center, an edge-like structure 
5' — 7' to the West from the center and complicated structure in the inner 1' — 2' region. 



Fig. 7 the relative deviations of the surface bright- 
ness from the azimuthally averaged surface bright- 
ness value (i.e. [Data — Model] /Model). The re- 
sulting image has been smoothed with a 2" Gaus- 
sian. Apart from the complicated structure in 
the innermost (~2') region probably dominated 
by the interaction of the thermal gas with the cen- 
tral AGN radio lobes (e.g. Bohringer et al., 1993, 
Churazov et al. 2000, Fabian et al. 2000, 2002) 
there arc two regions of prominent excess emission 
on large scales - one ~4'-7' west of the center and 
one > 5' east of the center. The same structures 
were also clearly seen in earlier ROSAT PSPC and 
HRI images (e.g. Schwarz et al. 1992). The possi- 
ble nature of these structures is further discussed 
in Section 4. 

3.4. Deprojection analysis 

The above radial dependencies of the gas pa- 
rameters are affected by projection effects and the 
azimuthal asymmetry of the cluster. The first 
problem can be solved using deprojection analy- 
sis. Various flavors of deprojection techniques are 
frequently used (e.g. Fabian et al. 1981, Kriss, 
Cioffi & Canizares 1983, McLaughlin 1999). In our 
analysis we assume spherical symmetry but make 
no specific assumption about the form of the un- 



derlying gravitational potential. We first calculate 
the surface brightness (in a given energy band) in 
a set of n a annuli and choose a set of n s spheri- 
cal shells. The radii of the annuli and shells need 
not be the same. The gas parameters are assumed 
to be uniform inside each shell. The outermost 
radii of the annuli and shells are chosen so that 
r(n s ) < r(n a ). The emissivity of the gas for radii 
larger than r(n s ) is assumed to decline with the 
radius as a power law with a given slope, namely 
£ = £$r ~ 4 . The normalization of this extra com- 
ponent So is an additional free parameter of the 
model. One can then write a simple expression 
which describes the contribution of all shells and 
outer regions from r(n s ) to infinity to the surface 
brightness S(j) in a given annulus j. 

S(j)= J2 A(i,j)E(i) + B(j)E , (1) 

i— l,n s 

where £(%) is the emissivity of a given shell. A 
simple analytical expression for A(i,j) and B(j) 
which are functions of the radii of corresponding 
shell and annulus can be easily written (e.g. Kriss, 
Cioffi and Canizares 1983, McLaughlin 1999). 

We then formulate a least squares problem - 
what set of emissivities in our set of shells will 
provide the best description of the surface bright - 
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ness recorded in our set of annuli: 
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where Oj is the statistical error associated with the 
surface brightness in a given annulus. As usual 
the differentiation of this relation with respect to 
S(j) and £ yields a system of linear equations, 
which can easily be inverted. The properties of 
the inverse matrix C(j,i) (in particular an error 
enhancement due to ill-conditioned problem) can 
be easily controlled by making the spherical shells 
broader. The emissivities of the shells are then 
expressed through the observed surface brightness 
distribution as an explicit linear combination: 



£[i)= ]T C(j,i)S(j) 



(3) 



We then accumulate a set of spectra (corrected 
for the background and vignetting) for each annu- 
lus. The emission spectrum associated with each 
shell can then be calculated as the same linear 
combination of the observed spectra. Since the 
whole procedure is linear the (statistical) errors 
can be propagated straightforwardly. The result- 
ing spectra for every shell are fit with standard 
models in XSPEC. The radial dependence of the 
electron density, temperature and heavy element 
abundance is shown in Fig. 8. As usual for any 
kind of deprojection analysis the errors on the de- 
projected values of temperature and abundance 
are significantly larger than for the projected spec- 
tra. However the biggest uncertainty in the de- 
projected parameters probably comes from the as- 
sumption of spherical symmetry. As is clear from 
Section 3.3 the cluster is significantly asymmet- 
ric. Furthermore when deriving the density in a 
given shell we assume that X-ray emitting gas uni- 
formly fills the volume of the shell. Nonetheless we 
present below two simple analytical approxima- 
tions to the density and temperature distributions 
which crudely characterize the radial behavior of 
these two parameters (shown in Fig. 8 with the 
solid lines): 
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Fig. 8. — Radial profiles of deprojected parame- 
ters. The parameters were obtained using a sin- 
gle temperature APEC model with fixed low en- 
ergy absorption fitted to deprojected spectra for 
a set of spherical shells centered at NGC 1275. 
The data points in the innermost region (~10 kpc) 
are affected by the presence of the bright compact 
source - the nucleus of NGC 1275. Analytical ap- 
proximations for the density and temperature are 
shown as solid lines (see text). 



where the second term is the density distribution 
on larger scales taken from Jones & Forman (1999) 
and shown in Fig. 8 with the dotted line. 
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(5) 



The analytical approximations shown in Fig. 8 
clearly fail to follow the data in the innermost 
region (within the central ~10 kpc). Given the 
XMM-Newton angular resolution, the major prob- 
lem here is the contribution from the bright com- 
pact source - the nucleus of NGC 1275. This 
makes Chandra much more suitable for studies of 
the innermost region. 

Schmidt et al. (2002), using Chandra data 
and excluding the contribution from the nucleus, 
showed that within about 50 kpc, the deprojected 
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temperature profile is essentially flat (see their 
Fig. 7) with a mean of approximately 3.1 keV in 
excellent agreement with the XMM-Newton tem- 
perature values from 10-50 kpc (away from the re- 
gion in which the nuclear source contributes to the 
XMM-Newton spectral data). The temperature 
profile derived from the projected Chandra spec- 
tra in finer radial bins (Fig. 2 from Schmidt et al. 
2002) does however show a slight increase in tem- 
perature in the innermost (<10 kpc) bin. It is pos- 
sible that projection effects (due to the presence of 
X-ray "holes" in the inner region) increase some- 
what the apparent gas temperature there. There- 
fore, equation 5, with its essentially constant tem- 
perature at r < 50 kpc is a reasonable approxi- 
mation. A comparison of the azimuthly averaged 
Chandra surface brightness profile in the inner re- 
gion with the XMM-Newton data also shows that 
a gas density profile with a flat core (eq. 4) is a 
reasonable approximation to the data. 

3.5. Projected Temperature Structure 

To calculate the projected gas temperature dis- 
tribution, we employed the method described in 
Churazov et al. (1996, 1999) and applied pre- 
viously to the ASCA data (e.g. Donnelly et al. 
1998, 2001). Briefly we use template spectra 
corresponding to emission from an optically thin 
plasma (convolved with the MOS energy response) 
with temperatures 3 and 9 keV with a given metal- 
licity and low energy photoelectric absorption and 
determine the best fit weights of these template 
spectra needed to describe the spectrum S(E) ob- 
served in a given pixel of the image, i.e. 

S(E) = A-M(T U E) + B-M(T 2 ,E), (6) 

where M(T,E) is the template spectrum for a 
given value of temperature T. The temperature 
is then calculated as a function of the relative 
weights of the template spectra. As is shown in 
Churazov et al. (1996) an expectation value of 
the temperature calculated this way is usually very 
close (within a few percent) to the temperature ob- 
tained by conventional spectral fitting under the 
assumption of a single temperature plasma with 
fixed metallicity and absorption. Compared to 
hardness ratios, typically used for the same pur- 
pose, this method is significantly more sensitive. 
The determination of weights is a linear proce- 
dure (i.e. fast) and further smoothing can be 



applied after the images containing weights are 
calculated. The resulting temperature map for 
Perseus is shown in Fig. 9. An adaptive smoothing 
has been applied to the map in such a way that 
each value of the temperature has been calculated 
using regions containing <~ 10 4 counts. Compar- 
ison with the results of direct spectral fitting of 
individual regions shows good agreement with the 
overall temperature structure. 

Compared with the previous ASCA tempera- 
ture maps (Arnaud et al. 1994, Furusho et al. 
2001) XMM-Newton data show the innermost re- 
gion in much greater detail, while the global asym- 
metry is broadly consistent with the lower resolu- 
tion ASCA maps, which cover a much larger area. 
The temperature structure of the very central re- 
gion (within 1' — 2') is consistent with the recent 
Chandra results (Fabian et al. 2000, Schmidt et 
al. 2002). 

On the spatial scales of 1-15 arcminutes, which 
are particularly well mapped by the XMM-Newton 
observations, the most remarkable feature of the 
temperature map is a horseshoe-shaped high tem- 
perature region to the West, enveloping the central 
cooler region. Thus the global east-west asymme- 
try of the cluster seen in the surface brightness 
distribution (Fig. 7) is also clearly present in the 
gas temperature distribution. 

If the cluster gas is in a state of a hydrostatic 
equilibrium then the gas density and temperature 
should remain constant along the equipotential 
surfaces set by the gravitating mass. In the exam- 
ple of a spherical or elliptical potential, all features 
in the observed surface brightness or temperature 
distribution should reflect the symmetry of the un- 
derlying potential. Shocks or strong sound waves 
(e.g. produced by a merger) may reveal themselves 
as distinct features in the observed images or tem- 
perature maps. The variations of density and tem- 
perature in shocks are obviously correlated leading 
to a positive correlation of the variations in surface 
brightness images and temperature maps. The ob- 
served variations (Fig. 7 and Fig. 9), on the con- 
trary, suggest that in Perseus the variations of the 
temperature and surface brightness are anticorre- 
lated. One can see that all the most prominent 
positive features in Fig. 7 correspond to regions of 
lower (relative to the mean value at the same dis- 
tance from the center) temperature. To show this 
more clearly in Fig. 10 we plot side by side the de- 
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Fig. 9. — Gas temperature map (30' x 30'). The temperature is derived from adaptively smoothed data with 
an effective number of counts within each smoothing window of ~ 10 4 . Contours show the surface brightness 
distribution. 



viations of the surface brightness and temperature 
from the azimuthally averaged values. The sur- 
face brightness image I and the temperature map 
T can be further combined in the form of two ad- 
ditional images P = T\fl and S = T/I 1 / 3 . Since 
the surface brightness scales as the square of the 
electron density (i.e. / cx n^) then P oc Tn e and 
S oc T/nl^ 3 . Therefore these two images P and S 
can be used to approximately characterize the az- 
imuthal variations of pressure and entropy 5 of the 
gas. These variations are also shown in Fig. 10. 
It is clear that the relative amplitude of the az- 
imuthal variations in pressure is smaller than that 
for any of the other quantities shown in Fig. 10 
and there is no clear east-west asymmetry in the 
pressure distribution. This argues against the at- 
tribution of the bulk of the observed substructure 
in the surface brightness distribution and temper- 
ature to strong shocks. Instead the observed fea- 
tures clearly resemble edges or "cold fronts" found 
recently in several clusters (e.g. Markevitch et al. 
2000, 2002, Vikhlinin et al. 2001). The possible 
origin of this substructure is discussed below. 



5 we use the term entropy for the quantity S = T/nJ , 
true specific entropy is determined by the logarithm of this 
quantity 



4. Discussion 

During the lifetime of any cluster, it experi- 
ences a few major mergers with another cluster 
of comparable mass and much more frequent ac- 
cretion of smaller clusters or groups. A chain of 
bright optical galaxies extending to the west from 
NGC 1275 suggests that at present the Perseus 
cluster is merging with a subcluster coming from 
that direction. The effect of a merger on the ap- 
pearance of the cluster obviously depends strongly 
on the mass ratios and gas fractions of merging 
clusters (see recent reviews by Evrard & Gioia 
2002, Sarazin 2002, Schindler 2002 and references 
therein) . Given that the deviations from spherical 
symmetry in the surface brightness and tempera- 
ture, although clearly present, are not very large 
then we probably have a case of a merger with 
a relatively small subcluster or group of galaxies. 
Furthermore the lack of large azimuthal pressure 
variations (Fig. 10) suggests that we observe the 
merger in a progressed stage when strong shocks 
have already left the central region of the clus- 
ter mapped with XMM-Newton. The numerical 
(N-body and fluid dynamics) simulations of the 
merger most relevant for the Perseus case are prob- 
ably those of Gomez et al. (2002). Their merger 7 
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Fig. 10. — Relative deviations from azimuthly av- 
eraged values of a) surface brightness I b) tem- 
perature T c) pressure p = T\fl and d) entropy 
s = T/I 1 / 3 . Contours mark the regions of ex- 
cess surface brightness (from the first image) . The 
color scale is the same for all three plots (from - 
0.5 to +0.5, see color bar in Fig. 7). Clearly the 
azimuthal variations in pressure are smaller than 
in any other quantity. 

(see Table 1 from Gomez et al., 2002) has a mass 
ratio 1:16 and small gas fraction for the subcluster 
so that the cooling flow is not disrupted during the 
merger. One can see some features in the gas den- 
sity and temperature distributions in their simu- 
lations (Fig. 5 and 6 in Gomez et al., 2002) which 
are qualitatively similar to the structures found in 
the Perseus data. 

In Perseus the most striking feature is a clear 
edge in the surface brightness which separates 
cooler gas in the very core from the hotter 
horseshoe-shaped region in the West, envelop- 
ing the core. One plausible explanation of these 
structures is a contact discontinuity separating 
main cluster gas from the gas of the infalling sub- 
cluster. The shape of the contact discontinuity 
depends primarily on the density profiles of col- 
liding clusters. For two identical clusters involved 
in a head-on merger the contact discontinuity is 
a plane perpendicular to the line connecting the 
centers of two clusters and located half way be- 



tween the centers. In the case of a merger of a rich 
cluster with smaller and less dense subcluster the 
contact discontinuity shifts towards the center of 
the smaller subcluster and bends due to stronger 
impact of the ram pressure onto the outer regions 
of the smaller subcluster, giving the infalling gas 
a conical shape (see e.g. Fig. 5 in Gomez et al. 
2002). The presence of a very dense core in the 
main cluster may reverse this trend and bend the 
contact discontinuity in the opposite direction so 
that it will envelop the very core of the main 
cluster. This behavior is less obvious in the simu- 
lations by Gomez et al. (2002) although flattening 
of the contact discontinuity is definitely seen. As- 
suming that merger 7 of Gomez et al. (2002) is a 
close analog of the Perseus merger, one would con- 
clude that the present state of the Perseus cluster 
corresponds to a time ~ 0.25 Gyr after the core 
crossing when the contact discontinuity is located 
within ~ 200 kpc of the center of the main cluster 
(upper right panel of their Fig. 4). After 0.5-0.75 
Gyr the dense and cold gas of the main cluster 
rebounds and shifts the bulk of the subcluster gas 
to larger distances from the center (middle right 
panel of their Fig. 4). 

An alternative explanation of the edge and 
horseshoe structures in Perseus assumes that the 
contact discontinuity is located farther away from 
the cluster core and the observed structures are 
solely composed from the disturbed gas of the 
main cluster. If the chain of bright galaxies 
marks the ongoing accretion direction of very 
small groups of galaxies, it is possible that (given 
the richness of Perseus) in each individual merger 
event the gas of the infalling subcluster is decel- 
erated far from the Perseus core and never comes 
close to the center. Such a situation resembles 
the picture of a piston which initially is moving 
through the gas but then stops. The gas in the 
core of the Perseus cluster would then see a com- 
pression wave (shock) followed by a rarefaction 
wave. The gas density and temperature increase 
during the compression phase, but then decrease 
as the gas reexpands. There should be only a 
modest increase in the final gas temperature as- 
sociated with the increase of the entropy in the 
shock. The velocity of the gas in the core of the 
main cluster also first increases and then falls back 
to zero, with the net effect that every volume el- 
ement of gas is displaced some distance from its 
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initial position. Since the density, temperature 
and pressure change significantly across the core 
(particularly in the direction perpendicular to the 
direction of the compression wave motion) the im- 
pact of the passing compression/rarefaction wave 
will differ in terms of the increase in temperature 
and the net displacement. If the characteristic 
amplitude of the wave propagating through the 
nonuniform medium scales as the square root of 
the density n e then the amplitude of the initial 
displacement at a distance r from the center X(r) 

scales as . / — \-r . Therefore the wave would in- 

duce a smaller net displacement in the very core 
region (where the density is the highest), than in 
the region farther from the core. Such a process 
can lead to characteristic "horseshoe" structures 
in the temperature and density distributions that 
are qualitatively similar to those observed in the 
Perseus cluster. However the relative amplitude of 
the displacement X(r)/r oc ^ n ^ r y 2 is a decreas- 
ing function of radius for density profiles less steep 
than r~ 2 . Therefore the relative amplitude of the 
displacement in the very core of the cluster will 
still be strong. The displaced gas in the core will 
then try to restore hydrostatic equilibrium and os- 
cillate in the potential. 

We illustrate the bending of the shock and the 
induced oscillatory motion of the gas in the core 
with a simple numerical simulation of a plane com- 
pression wave passing through a cluster gas. The 
simulations were made using the ZEUS-2D code 
(Stone & Norman 1992a, b). In our simulations 
we employed an ideal gas equation of state with 
7 = 5/3. Equations 4 and 5 were used for the ini- 
tial density and temperature distributions. The 
gravitational potential was calculated from the 
same density and temperature distributions and 
assumed to be static. The computational domain 
spans a 6 x 4 Mpc region and is covered by 600 x 
400 grid points. At some distance from the cluster, 
we modify the gas velocity, density and tempera- 
ture in a slab 200 kpc wide, to mimic the situ- 
ation when a piston is first driven into a cluster 
gas with a Mach number of ~2 and then stops. 
The results shown below were obtained from the 
simulations on a Cartesian grid, while simulations 
on the cylindrical grid produce qualitatively very 
similar results. 

Snapshots of the gas density distribution at sev- 




Fig. 11. — Snapshots of the gas density at 0, 0.3, 
0.7, 1.6 Gyr from the beginning of the simulations. 
The box shown is 4x6 Mpc. The density, temper- 
ature and velocity of the gas in a narrow slab on 
the right side of the box were set to mimic the 
impact of a piston which is first driven into the 
gas at the Mach number of 2 and then suddenly 
stopped. 

eral successive moments of time (0, 0.3, 0.7, 1.6 
Gyr from the beginning of the simulations) are 
shown in Fig. 11. Initially the plane shock bends, 
when passing through the core, because of the 
higher pressure and lower temperature in the cen- 
tral region. Shortly after the passage of the wave 
through the core, the gas distribution is clearly 
elongated in the direction of motion, with the 
very core regions experiencing smaller displace- 
ments than the outer layers. The gas then begins 
"sloshing" (this term was introduced in Marke- 
vitch, Vikhlinin & Mazzotta 2001, see also Marke- 
vitch, Vikhlinin & Forman 2002) in the potential 
well, which in our case is assumed to be static. 
To show the structure of the density distribution 
more clearly, we plot the relative deviation of the 
gas density from the undisturbed value in Fig. 12. 
The first frame corresponds to the time when the 
shock front is passing through the core. Although 
originally the gas is swept in the direction of the 
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Fig. 12. — Relative deviation of the gas density from the initial value (blue color marks the regions with 
the density below the initial value while yellow and red correspond to overdense regions) . The box shown is 
approximately 2 Mpc on a side. The snapshots are taken at 0.3, 0.7, 2.3, 3.9 Gyr from the beginning of the 
simulation. 



wave propagation, in ~0.5 Gyr it rebounds and 
produces excess emission on the opposite side of 
the core. The relevant time scales for the gas 
sloshing are probably set by the acoustic cutoff 
frequency w a and the buoyancy (Brunt- Vaisala) 
frequency N, which can be written as 



10" 



N = fl K A 



c m d InP 
idr 



/ 1 d Ins 
7 d In r 



(7) 



(8) 



where s is the gas entropy, r is the radius, and Q,k 
is the Keplerian frequency at a given radius. The 
dependence of the acoustic cutoff and the buoy- 
ancy frequencies on radius (using the density and 
temperature profiles according to equations 4 and 
5) are shown in Fig. 13. 

While strictly speaking these frequencies are de- 
termined by the local properties of the gas (and 
we are interested in global oscillations), it is likely 
that they set an approximately correct timescale 
for the gas sloshing. Since these frequencies vary 
with radius, the gas will slosh with different peri- 
ods, with the possibility of forming several asym- 
metric structures at the same time on either side 
of the core. This is illustrated in Fig. 12 where the 
deviation of the gas density from the undisturbed 
value is shown at several times. The structures 
seen in the density distribution resemble the cold 
edges or fronts observed by Chandra in many clus- 
ters (e.g. Markevitch et al. 2002). The edges in 
this simulation are due to sloshing of the cluster 
core gas and are of course transient events, which 
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Fig. 13. — Buoyancy N (solid line) and acoustic 
cutoff w a (dotted line) frequencies as a function of 
radius. 



however can persist for a few Gyr, although their 
shapes evolve continuously. As usual for oscilla- 
tory motion, the gas spends most of the time near 
one of the points of maximum displacement - i.e. 
in the configuration favorable for the detection of 
the disturbed gas as an "edge" . 

From Fig. 13, the shortest period of sloshing is 
on the order of a Gyr. This means that it takes 
about a 0.5 Gyr for the gas displaced by the ini- 
tial shock to rebound to the opposite side of the 
center of the potential and form an edge there. If 
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this mechanism is responsible for the formation 
of the edge in the Perseus cluster then we can 
roughly estimate that the initial shock has passed 
through the core some 0.5 Gyr before. On such 
time scales the primary wave has enough time to 
propagate outside the region mapped with XMM- 
Newton and the disturbed dense gas in the center 
has enough time to rebound. 

In this interpretation the increase of gas entropy 
on the western side of the cluster (see Fig. 10) can 
be attributed to the initial shock or sequence of 
shocks. If the gas in that region was originally on 
the same adiabat as now observed at the eastern 
side of the cluster then in a single shock inter- 
pretation the increase of entropy by ~50% (esti- 
mated from Fig. 10) would require a shock with a 
Mach number of ~2.6. Assuming that the grav- 
itating mass profile of the Perseus cluster follows 
the Navarro-Frenk- White profile (Navarro, Frenk 
& White, 1997) with the parameters obtained by 
Ettori, De Grandi and Molendi (2002) we can es- 
timate that the free fall velocity reaches a Mach 
number of 2.6 (for 6.5 keV gas) at a distance 
of ^300 kpc from the cluster center. This is of 
course an upper limit on the velocity of the con- 
tact discontinuity, attainable only if the infalling 
gas is sufficiently dense. Depending on the ge- 
ometry of the contact discontinuity the forward 
shock will have a velocity similar to the velocity 
of the discontinuity (e.g. bow shock) or slightly 
higher (for a plane parallel shock). Therefore the 
observed asymmetry in the entropy distribution 
is only marginally consistent with the assump- 
tion that a single shock was responsible for the 
entropy increase. Multiple small accretion events 
from roughly the same direction, which can be ex- 
pected if there is indeed a large scale filament in 
the direction set by the chain of bright galaxies 
could more easily provide the observed rise of the 
entropy. 

The motion of the cD galaxy relative to the 
larger scale potential of the cluster has long been 
suspected to play a role in clusters with cool and 
dense cores (e.g. David et al. 1994, Fabian et 
al. 2001, Markevitch et al. 2001). Since the dis- 
tinction between the cD potential and cluster po- 
tential may be somewhat artificial, one can think 
instead in terms of perturbations of the potential 
in the cluster core. The gas can then start sloshing 
in the disturbed potential (e.g. Markevitch et al. 



2001, 2002). We argue above that pure hydrody- 
namics may cause similar sloshing. What matters 
is the displacement of the gas relative to hydro- 
static equilibrium, which can be achieved either 
through disturbances in the gas or the gravitat- 
ing mass. In real mergers, both mechanisms are 
probably operating. 

Markevitch et al. (2001, 2002) suggested that 
dissipation of the sloshing kinetic energy may be 
sufficient to offset radiative gas cooling. Assum- 
ing that the asymmetry in the surface brightness 
distribution in the inner 200 kpc region around 
NGC 1275 directly characterizes the amount of 
sloshing energy, one can estimate this energy as 
10-20% of the thermal energy of the gas. At a 
radius of 100 kpc, the cooling time is comparable 
to the sloshing turn around time (P = 2-K/w a or 
P = 2n/N) and at 200 kpc it is several times 
longer than the sloshing time. Assuming that 
this energy is dissipated over several sloshing turn 
around times (and neglecting the fact that part 
of the energy would be emitted as acoustic waves 
which escape the region) , we conclude that the dis- 
sipated energy is not sufficient to offset cooling. 
However, the magnitudes (for cooling and dissi- 
pation) are comparable and one can imagine sce- 
narios where the dissipation of the sloshing kinetic 
energy plays a more important role. In particular 
our estimate that the sloshing energy constitutes 
10-20% of the thermal energy may be an underes- 
timate. 

4.1. The fate of the stripped gas 

During the merger process the subcluster gas 
is stripped from the dark matter potential of 
the infalling subcluster and mixes with the main 
cluster gas. As was pointed out by Fabian & 
Daines (1991), the infalling gas does not neces- 
sarily pass through a very strong shock during the 
merger process and therefore the gas may remain 
on approximately the same adiabat as before the 
merger. This means that even stripped and de- 
celerated lumps of subcluster gas may continue 
to move towards the center of the main cluster 
if their entropy is low enough. Even neglecting 
mixing with the main cluster gas, this motion can 
continue only until the lump reaches the radius 
in the main cluster where the main cluster gas 
has the same entropy as the lump. In Fig. 14, 
we plot the entropy (defined as T e /n^ 3 ) profile 
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Fig. 14. — Radial gas entropy profile calculated as 
T e /n^ 3 using the deprojected density and tem- 
perature profiles. 



of the gas in the Perseus cluster. One can com- 
pare the values of the entropy with the character- 
istic values of the entropy found in clusters and 
groups (e.g. Ponman, Cannon & Navarro 1999; 
Lloyd-Davies, Ponman, & Cannon 2000) which ex- 
ceed ~100 keV cm 2 even for poor groups. Thus, 
even if the infalling gas does not pass through a 
strong shock, which would further increase its en- 
tropy, it will find gas with comparable entropy at 
a distance at least several hundred kpc from the 
Perseus core. The infalling gas cannot easily pen- 
etrate below this radius. Therefore the core of the 
Perseus cluster is well "protected" against pene- 
tration by infalling gas (unless the infalling sub- 
cluster itself contains a cool core with low entropy 
gas). 

4.2. Mergers in the absence of a dense 
(cool) core 

The picture discussed above assumes that the 
main cluster contains a very dense (cool) core, 
while the infalling subcluster does not. In real- 
ity, various combinations of these parameters are 
possible, leading to a wealth of different appear- 



ances for mergers. In particular, if the main clus- 
ter lacks a dense (cool) core and the entropy of the 
infalling gas is significantly lower than that of the 
main cluster, the infalling gas may penetrate deep 
into the cluster core. The infalling cooler gas will 
be progressively stripped (Fabian & Daines 1991) 
from the dark matter halo and partly mixed with 
the main cluster gas. Even when stripped, the gas 
remains overdense compared to the main cluster 
gas and it will keep moving (although more slowly) 
towards the bottom of the potential well. This 
would result in a "cool wedge" penetrating deep 
inside (or even crossing) the main cluster. This 
kind of structure is probably observed in A1367 
(Forman et al. 2003) and in Coma. Both clusters 
have very high values for the gas entropy in their 
cores, so it is relatively easy for infalling gas to 
penetrate the core. Thus the cores in such clusters 
may appear much more disturbed than the cores 
of the "cooling flow" clusters, even if the infalling 
subclusters have similar parameters. 

Thus if the main cluster does not have a dense 
and cool core, the gas from smaller infalling sub- 
clusters may penetrate deep into the cluster and 
mix with the main cluster gas. One can there- 
fore expect the abundance of heavy elements to 
be more or less uniform throughout the cluster. 
However when the main cluster has a dense cool 
core, the infalling gas is stopped before reaching 
the core and only small lumps from any cool cores 
in the infalling groups themselves arc "allowed" 
to penetrate deeper. These infalling cool lumps 
are generally centered on a bright galaxy and are 
presumably also enriched in metals compared to 
the bulk of the subcluster gas which is stopped at 
larger distances from the center of the main clus- 
ter. Therefore this process may serve to maintain 
(or even enhance) the radial abundance gradients 
observed in many clusters with cool cores. 

5. Conclusions 

We present the first results of the 50 kscc XMM- 
Newton observations of the Perseus cluster. The 
surface brightness and the gas temperature maps 
show a wealth of substructure on all spatial scales. 

The inner region (1' — 2') is clearly disturbed by 
the interaction of the outflow of relativistic plasma 
from NGC 1275 with the thermal gas. On larger 
scales the horseshoe structure of the gas tempera- 
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ture and the surface brightness structure suggest 
an ongoing merger along the direction defined by a 
chain of bright (optical) galaxies. The gas from the 
infalling subcluster may be stripped from the dark 
matter halo at a large distance from the core of the 
Perseus cluster, with only compression/rarefaction 
waves passing through the core of the main clus- 
ter. These waves may induce "sloshing" in the gas 
in the cluster core. As a result an "edge" (or even 
several edges) in the surface brightness and in the 
gas temperature distribution may appear. 

Similar structures may be present in other clus- 
ters with peaked (cool) cores which are presently 
accreting smaller subclusters (e.g. Centaurus clus- 
ter). For clusters without a cool dense core, the 
situation may be different since the subcluster gas 
is cool and can penetrate deep inside the core of 
the main cluster and form a cool wedge of partially 
stripped and mixed gas (e.g. Coma, A1367). 
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